NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



((IASA-Tn-d2169) INIEFPCLAIING FOB IHE 
LOCATIOB OF BEtiOTE SENSOfi tATA (BASA) 23 p 
UC A02/BF A01 CSCL 05B 


B82-13469 


Uncias 
G3/43 04371 


IVI/>SA 

Technical Memorandum 82169 


Interpolating for the Location of 
Remote Sensor Data 


E. F. Puccinelli and J. Kornfield 




JULY 1981 


National Aeronautics and 
Space Administration 

Goddard Space Right Canter 

Greenbelt. Maryland 20771 




INTERPOLATING FOR THE LOCATION OF 
REMOTE SENSOR DATA 


by 

Edward F. Pucdnelli 
NASA/Goddard Space Flight Center 
Code 931 

Greenbelt, Maryland 20771 
and 

Jack Komfield 
National Science Foundation 
1800 G. St. 

Washington, D.C. 20SS0 


July 1981 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 



INTERPOLATING FOR THE LOCATION OF 
REMOTE SENSOR DATA 

Edward F. Puccinelli 
NASA/Goddard Space Flight Center 

and 

Jack Komfleld 
National Science Foundation 

ABSTRACT 

An interpoladon algorithm is presented as a practical alternative to common inter* 
polation and appMximation methods when applie^ to the problem of determining the 
location of remote !«nsor data. This algorithm is based upon knowledge of the geometry 
of the problem and is shown to be inherently more accurate than common interpolation 
schemes which may be applied to all types of data. A practical location problem is used 
to demonstrate its accuracy and computational cost. 
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INTERPOLATING FOR THE LOCATION OF 
REMOTE SENSOR DATA 


INTRODUCTION 

A common problem encountered in the analysis of remote sensor data is the determination 
of the location of large volumes of data. The number of computations for the location of any 
single field of view is in itself not great. When this number is multipli^xl by the high sampling 
rate common to many of today’s scanners, the cost in terms of time and money may become 
intolerable. This fact has lead analysts to turn to the use of interpolatiun and approximation 
methods to aide in determining the location of fields of view at a reasonable cost. Usually, ac* 
curacy of location is sacriHced for speed of computation. 

First, this paper will investigate the possible reasons for location inaccuracy which arise when 
using interpolation or approximation schemes. The results of the discussion will be used to intro- 
duce an interpolation algorithm which serves as a practical alternative to the more traditional 
methods of spline interpolation, polynomial interpolation, least squares polynomial approximation, 
etc. This new algorithm will then be compared to spline interpolation in an actual application. 

BACKGROUNDS, ASSUMPTIONS AND DEFINITIONS 

A simple example serves to introduce the interpolation problem. A remote sensor will be 
assumed to be on board a spacecraft orbiting the earth. The discussion which follows is not de- 
pendent on this assumption and the conclusions are equally valid for sensors on board aircraft or 
for bodies other than the earth. 

Other assumptions that will be made for the sake of clarity of the discussion are associated 
with the sensor and the spacecraft orientation. The spacecraft axes will be considered to have 0** 
pitch, 0° yaw and 0° roil attitude when the yaw axis is normal to the ellipsoid representii^ the 
earth surface and the pitch axis is coincident with the cross product of the normal to the surface 
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and the spacecraft velocity vector. The sensor axes will be assumed coincident with the space* 
craft axes. Figure 1 illustrates the geometry under these assumptions. 



It is also assumed that the sensor rotates in the pitch-yaw plane so that in one rotation the 
sensor returns data for n fields of view. These n fields of view will be defmed as being one scan 
of the sensor. The fust field of view of the scan will be defined as occuring at time t{ and the 
last at time t„. ^ assumption is made about the angular rate of motion of the sensor about the 
roll axis or about the sampling rate. 

Again, these assumptions made for the sake of clarity do not restrict the generalization of 
the following discussion to other types of sensors. Figure 2 illustrates the second set of 
assumptions. 

To lay a framework for discussing interpolation errors it is necessary to establish the infor- 
mation needed in order to locate a single field of view. The following is a list of the data 




Figure 2. A Single Scan as Viewed by a Plane 
Rotating Sensor on Board a Stationary Spacecraft. 


necessary for tliis task. (All information, except when otherwise spedfled is measured relative to 
an earth fixed cartesian coordinate system whose origin is coincident with the center of the earth.) 

1. An equation describing the earth surface. C^ail this E. 

2. The position of the spacecraft as a vector function of time. Call this'stt). 

3. The velocity of the spacecraft as a vector function of time. (Tall this't^t). 

4. The aritude of the spacecraft as a vector function of time. Call thisT(t) where*a(t) » 

(pitch (t). roll (t), yaw (t))^. 

5. A vector function iT(t) which describes the sensor pointing direction at time t relative to 
the sensor axes. (Recall that the sensor axes and the spacecraft axes are assumed coinci* 
dent for the sake of simplicity. Were they not. their relationship would be a sixth piece 
of information.) 
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Reference 1 describes in detail a method of determining a vector function f of these five 
parameters which yields the earth fixed locationTof sensor data. 

■ctt)-T(E,ttt),Tl[t).-r(t),tr(t)) (1) 

If the four functions of time (ttt),T(t),'3Xt), ^t)) and equation E were known exactly at 
any time t, then the vector *^t) would be the exact location of the data sensed at t. If from 
scan start t^ to scan end tg the input functions were continuous thenT would yield a continuous 
trace of the scan on the earth surface E. 

Evaluation ofTat discrete times tj, t 2 , . . . , tg when the sensor is recording data will pro- 
duce discrete points which are the locations of the individual fields of view in a single scan. This 
set of locations {'T(tj),'e(t 2 ), . . . ,’^(tg)} is defined as a scan trace. 

Of course, these five input parameters are hardly ever known accurately which in itself leads 
to error in the determination of any single location. An analysis of this type of error is beyond 
the scope of this paper. The discussion here is confined to developing an interpolation method 
which permits rapid determination cf the location of the scan trace. 

THE INTERPOLATION PROBLEM 

In practice when locating remote sensor data the usual procedure is to determine several lo- 
cations in a scan trace as accurately as possible. This subset of locations will be referred to as 
anchor points.* An interpolation or approximation algorithm is then applied to the anchor points 
in an effort to determine the location of the remainder of the sci.n trace. Ignoring the errors in 
location of any individual field of view due to imprecise or inaccurate measurements of the five 
necessary inputs, the location problem can be expressed as follows: 


*The anchor pointi may be expreoed in eiUier cartesian coordinates or latitudes and longitudes. Usually, it is computationally advan- 
tageous to work with latitudes and longitudes, however, for the ake of being able to envision the problem in three dimensions the 
discussion will refer to catted coordinates. A pnctical example is presented later for both systems. 
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Given a set of anchor points determine a vector fiinction f such that 

ff(t)-T(t)l2<« (2) 

for all t« {ti, tj, . . . ,tn}. 

The vector norm I T 1} is defined as the Eudidean length of the vectorT 

Note that expressing the problem in this fashion allows for the use of either interpolation or ap* 

proximation methods for Ending'S! 

SOURCES OF INTERPOLATION ERRORS 

There are two main reasons why lf(t) - g(t)l 2 may exceed 6 between anchor points'^ • 
"etti) and'^ ='f(tj). 

First, *g1[t) may iu>t coincide with the earth surface. An exaggerated but simple example will 
demonstrate this point. Suppose T(t) is the straight line through two anchor points'^ and*^ and 
that the sensor scan plane does not contain the earth center but is actually sensing data near the 
earth horizon. Figure 3 shows how the difference between f andX might be larger than a sped- 
Eed tolerance. 

It may be argued that if approximation or interpolation is performed on latitude and longi- 
tude coordinates then the resulting locations are implidtly on the earth surface. However, as 
Figure 3 shows, if the scan plane is not near the earth center there may still be a significant error 
in approximating f. 

The second prindple source of error is dependent on the behavior of the four functions 
's(t),'v(t),T(t),^(t). Respectively, these describe the position, velodty and attitude of the space- 
craft and the motion of the seiuor relative to its axes. Consider the following idealized example. 

1 . Let E describe the earth surface exactly. 

2. Let (tj, tj) be the time interval in one scan during which locations are sought. 

3. Assume that the spacecraft does not move during (t^, tj) and is exactly over the equator 
and the 0° longitude line. 
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SENSOR 


4 Assume that the velocity vector orients the spacecraft roll axis so that it points due 
north. 

S. Let the scanner pointing vector rotate at a constant rate about the roll axis in the pitch* 
yaw plane. 


Figure 4 shows how the scan trace will coincide with the equator under these conditions if 
the attitude were perfectly zero in all axes over (tj, tj). This figure also shows what the scan 
trace would be if roll and yaw were constantly zero but pitch varied sinusoidaly as described by 


pitch (t) 



If only two or three anchor points are used in the interval (t|, tj) the scan trace for constant 
attitude is well approximated but there is little hope for a good approximation in the case of 
varying pitch attitude. 
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EARTH 

CENTER 



When the idealized restrictions are removed and the spacecraft position and velocity are 
allowed to vary with time as well as all three attitude parameters it is not difficult to envision 
very complicatec( scan traces on the earth’s surface. 

The trick is to select a sufficient numbet of anchor points in the interval (t|, tj) to allow for 
accurate interpolation or approximation. In order to get an idea of what is a sufficient number 
the analyst must know the expected behavior of the functions'stt).'^t),T(t) and Tv (t). 

The usual approach to solving the basic problem described by inequality (2) is for the analyst 
to solve explicitly for m anchor points in a scan of n fields of view and use them to produce an 
interpolation or approximation functionT Common methods include polynomial and spline in* 
terpolation as well as least square approximations. Linear interpolation is in itself computarion- 
ally the least expensive of these methods. However, accuracy requirements may necessitate the 
computation of so many anchor points that it may be less costly to use another method which 
requires fewer anchor points to produce the sam<* accuracy. 
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An important point to note is that all of these methods attempt to approximate the scan 
trace functionT without using any information about the component functions ofT 

This fact leads to the concept of developing an interpolation algorithm which approximates 
the component functions rather thanTitself. Using these approximate values,Tis then evaluated 
exactly. 

AN ALTERNATIVE TO COMMON INTERPOLATION AND APPROXIMATION METHODS 

This new algorithm is described as follows: 

1. Within the time interval of a single scan (t^, 1^) select a time interval (t{, tj) for which 

rT(E,t(t),^t),T(t). W(t)) < 6 (3) 

where^,'^ anda^ are constant approximations of the spacecraft position, velocity and 
attitude and1T(t) is an approximation in an earth fixed cartesian coordinate system of 
y/{t) the sensor pointing direction function. TheM approximations are described in the 
following steps. Note that determining (t|, tj) ;nay well be impossible to do precisely 
but is in fact equivalent to selecting anchor points for any of the previously mentioned 
interpolation and approximation methods. 

2. The approximation to the position, velocity and attitude are 

Sc “ s(th) 

Vc » v(th) 

3c » a(tjj) 

where 

th • (ti + tj)/2 

3. Next, the vector function ~w(t) must be approximated. First determine the anchor points 

and*^ by evaluating ffE.'fttil.'^tt^l.’altjL'^fti)) for d « i, j. Then form the two 
unit vectors 

ki »(td -1h)/*r^ 
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These two unit vectors nuy be pictured at orifinttiiii froni the constant spacecraft posi- 
tion and pointing toward the two anchor points. Using these two vecton and knowl- 
edge of how the sensor moves with respect to the spaceitaft axes the functk>n1T(t) is 
required to approximate the motion and satisfy 

T(td) -lEJi for d • i, j 
nr(t)lj ■ 1 for t e (t|, tj) 

Further definition oflT(t) is dependent upon the motion of the sensor of interest. .\n 
example of the complete definition of k(t) for a sensor which rotates about the space- 
craft roU axis is given in the next section. 

4. Ihe final step of this algorithm is to evaluate 

ttt) »T(E,t;,‘rc.Tc,T(t)) for t e (tj, tj). 

This is done by finding the intersection of the line in space which approximates the sen- 
sor pointing direction with the earth surface. The line is given by 

Uu) uWt) 

and the earth surface by 


( a cos Ui cos U 2 
a sin U( cos U 2 
c sin U 2 

Equating Uu) to E(U}, U 2 ). the parameter u may be solved for by multiplying the third 
equation by a/c, then squaring b«nh sides of the equation, and finally adding all three 
equations together. This eliminates the parameters u^ and U 2 and gives a quadratic equa- 
tion in u. 

Au^ 2Bu C ■ 0 


where 


q ■ a^/c* 

A - k*(t) * k*(t) qkj(t) 
® * *xh*^(*) ”*■ 
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(4) 

(5) 

( 6 ) 


It should be noted that soiutioiM for u represent the dittance from the spacecraft to the 
intersection(s) of L(u) with the earth’s surface. Therefore the smaUCT value is the one of 
interest giving 

- AC 

u ■ ' 

A 

The situation is shown in Figure 5. The locationTis then evaluated as 

+iS(t). 



In sununary. the proposed algorithm involves 

a. selection of anchor points in a scan 

b. assumption of a single spacecraft position velocity and attitude between any two anchor 
points. 

c. the formation of a function k(t) which approximates the sensor motion between the two 
anchor points. 

d. determination of the intercept of k(t) with the earth surface. 
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The basic concept is that the exact function~?tB,'?(t),T(t),T(t),'w(t)) is being appn'ximated 
by making approximations of its component functions. 

Incidentally, the algorithm may also be useful over time intervals (t{, tj) where the attitude 
is not approximately constant. For instance, for a scaimer rotating in a plane about the space- 
craft roll axis, the roll attitude rate should be approximately constant. The pitch and yaw atti- 
tudes should be approximately constant. The motion of the spacecraft about the roll axis will 
then be represented by the scanner pointing vector 1T(t). 

A COMPARISON OF SPLINE INTERPOLATION AND THE PROPOSED ALGORITHM 

In order to evaluate the proposed algorithm, a count of the operations needed to locate one 
scan trace will be made. A similai' count will be made for the spline interpolation technique. 

After a comparison of efficiency is made the results of a test to compare accuracy will be 
described. 

L^t m, be the number of anchor points needed in order for cubic spline interpolation on 
latitude and longitudes to produce locations to an accuracy of 5. Let n, be the remaining number 
of fields of view in the scan to be located by the spline interpolation, (e.g. m, + n, => total num- 
ber of fields of view to be located in one scan). Let lUp and np be similar numbers for the pro- 
posed algorithm. Table 1 shows the number of operations necessary to determine the location of 
any single anchor point when evaluating the function f as described in reference 1. 

Table 2 shows the number of operations involved in each step of the procedure when locat- 
ing a scan trace of cubic spline interpolation of latuude and longitude coordinates. 

a. The longitude values must be checked and adjusted as necessary to avoid discontinuity 
problems in steps 2 and 3. These operations are not included in the count. 

b. The number of operations was determined from the code of reference 2, specifically, 
subroutines ICSICV and ICSEVU. 
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The number of operations required by the proposed algorithm will be dependent upon how 
the function l^t) is formulated to approximate the motion of the sensor of interest. In order to 
present a comparison with spline interpolation we will proceed assuming the case of a sensor 
which rotates in a plane about the ipacecraft roll axis. The formulation ofT(t) for this case is 
presented as follows: 

1 . Form the unit vectors Iq and kj as described in the previous section. 

2. Compute the angle a between the two vectors. 

a » oor^ (iq • Iq) 

3. Let 



tj-ti 


4. Then 

T(t) » sin (a - 0)1q + sin /siq 

is a function which approximates '^t). Note that it is not of unit length as formally re- 
quired, but since its maximum length is less than 2 there is no ambiguity in solving for 
the intersection of L(u) and E. The function '0t) is not normalized in order to save 
computations. 

The number of operations for the proposed algorithm based upon the given formulation of 
l^t) is given in Table 3. 

In Table 3, operations which would only need to have been computed once, no matter how 
many scans were being approximated, were not counted. Other methods of formulating 1^t) 
might be found for a specific problem which involve fewer operations. 

Tables 2 and 3 were set up to produce locations in latitude and longitude coordinates. This 
is advantageous to spline interpolation when comparing the totals. Had it been required to pro- 
duce locations in cartesian coordinates then in Table 2 the counts irt step 1 would be slightly re- 
duced while the counts in steps 2 and 3 would be multiplied by a factor of 3/2 to reflect 
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interpolation on three coordinates as opposed to two. Table 3 would be altered by deleting step 
4. The comparison is shown in Table 4. Note that constant terms have been purposely dropped. 

Table 4 shows that it is not obvious which method is mon,' efndent without further informa* 
tion about actual values of m,, n,, mp, np, and the time each type of operation takes on a spe* 
dfic computer. In practice, one must determine the values m, and nip for a specific problem and 
then refer to the table to establish which method is faster. 

In order to compare the accuracy of the two methods, an exact computer model was de* 
veloped for the scan trace of the Temperature Humidity Infrared Radiometer (THIR) on board 
the Nimbus series spacecraft. The parameters of the Nimbus -6 sun synchronous orbit used by 
the model were 

a. 99.15° inclination 

b. 7333.16km semi-n)^or axis 

c. .001 eccentricity 

d. 104.16 min period 

The modeled characteristics of the sensor are identical to those previously assumed. In this 
case, the sensor rotates at 48 rpm in a plane about the spacecraft roll axis at a linear rate. Each 
field of view is sampled for .22ms with a delay of .98ms between samples (reference 4) so that 
an entire scan from earth edge is approximately comprised of 343 samples. The spacecraft axes 
were oriented as assumed previously so that for 0 ° yaw, pitch and roll, the yaw axis is normal to 
the earth ellipsoid and the pitch axis is normal to the velocity and yaw vectors. The earth ellip- 
soid was assumed to have an equatorial radius of 63'78.144km and a polar radius of 6356.759km. 
For the time period of one scan the attitude was held at zero in all axes. The model then eval* 
uatedTat 343 times. The solutions . . •‘^43 were considered to be exact. 

The location accuracy requirements were chosen to be 5 * 3.6 km. The location volume re- 
quirements were for 93 fields to be located in each scan such that they were approximately evenly 
spaced in distance on the earth surface. 
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A computer program was written which usedTii'f|72t'9)43 as anchor points. The cubic 
spline through those points was evaluated at the remaining 340 locations and compared to the 
exact solutions. The number of anchor points was then incremented by five and selected so that 
they were evenly spaced along the scan on the earth surface. Interpolated locations were again 
compared to exact locations. This procedure was repeated incrementing the number of anchor 
points by five each time until the maximum difference between the exact solution and the inte^ 
polated solution at any field of view was less than 3.5 km. The result was that 23 anchor points 
produced an error of 6km while 28 anchor points reduced the error to 3 km. The error measure- 
ment was made first for spline interpolation of latitudes and longitudes and then for cartesian 
coordinates. No significant difference in error measurement was formed. 

Next the proposed method was applied. It was found that only two anchor points e^ and 
^343 save a maximum error of .5 km! 

Applying m, » 28, mp » 2, n, » 65 and np » 91 to Table 4 produces the comparison shown 
in Table 5. 

The proposed method is more efficient when interpolating on latitude-longitude and cartesian 
coordinates. In both cases the maximum error was less by a factor of 6 when the proposed 
method was used. 

It should be noted that if all 343 fields of view had to be located rather than 93 evenly 
spaced in distance, then Table 4 would have shown spline interpolation to be faster. 

CONCLUSIONS 

The proposed algorithm for location of remote sensor data is a practical alternative to com- 
mon interpolation and extrapolation techniques. If applied to the same anchor points used for 
one of the traditional methods it can produce more accurate locations. An important feature is 
that it requires the determination of fewer anchor points. Its efficiency is dependent upon tht 
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number of anchor points required to produce a desired accuracy and the munber of locations 
sought and may be easily evaluated using Table 4. 

In practice, for situations where extremely large volumes of locations are desired a hybrid 
approach may prove advantageous. If one considers linear interpolation to be the fastest inter* 
polation possible then it would be logical to determine the anchor points necessary in a scan to 
produce locations to a specified accuracy using this method. However, instead of determining 
these anchor points directly by evaluating f(E,'STt),'^t),'?(t),'w(t)) use the proposed scheme 
or perhaps one of the traditional methods to determine the anchor points needed for linear 
interpolation. 



Table 1 

Operations in Locating One Andior Point 


Solution Coordinates 

Mults. 

Adds 

V— 

Trig. 

Trig.-1 

Lat-Lon 

95 

43 

5 

14 ■ 

3 

Cartesian 

90 

42 

4 

14 

1 


Table 2 

Number of Operations When Using Spline Interpolation 


Step 

Mults. 

Adds 


Trig. 

Trig.** 

1. Locate m, anchor pc>ints in 
lat-lon coordinates*. 

95 m, 

43 m, 

5 m, 

14 m. 

3 m. 

2. Determine the coefficients 
of the cubic splines for lat* 
itude and longitude between 
each pair of anchor points**. 

32 m, 
-38 

38 m, 
-48 

0 

0 


3. Evaluate the spline functions 
at n, times for latitude and 
longitude values**. 

6n, 

8 n, 

0 

0 

0 

Totals 

6 n, + 

127 m,- 38 

8n,+ 

81 m,-48 

5 m, 

14 m, 

3m, 


a. The longinide values must be checked and adjusted u necessuy to avoid discontinuity problenis in steps 2 and 3. Thesa opetations 
are not included in the count. 

b. The number of operattons was determined from the coda of reference 2. Specifkally subroutines ICSICV and ICSEVU. 
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Table 3 

Number of Operations When Using the Proposed Algorithm 


Step 

Mults. 

Adds 

xT" 

Trig. 

Trig.-! 

1. Locate nip anchor points in 
cartesian coordinates 

90 mp 

42 mp 

14 mp 

14 mp 

nip 

2. Compute the nip - 1 angles 
a 

3 (nip - 1) 

2(mp-l) 



Hi 

3. Form the nip vectors k 

4mp 

5 nip 

mp 



4. Form np values of 0 

■n 

"P 




5. Evaluate1T(t)*sin(a-|3) 
1q sin (^)Tcj np times 

6np 

4np 


2np 


6. Find the intersection of L(u) 
and E a total of n. times. 
(See equations 4. 3, 6) 

8 Up 

+ 4(mp-l) 

4np 

+ 3 (ntp - 1) 




7. EvaluaLe the location*T* 
-t* u0t) a total of np 
times 

3np 

3 np 




8. Convert the np locations to 
latitude and longitude values 

5np + 
5 mp 

«p +nip 

Hp +mp 


2 itp 
+ 2mp 

Totals 

24 np 
+ 106 ra- 
-7 

13np 
+ 53 m- 
-5 

"P 

+ 6mp 

2np 

14 mp 

2np 

+ 4mp 
- 1 
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Table 4 

Compariaon of the Total Operations 


Coordinate System 
Desired for 
Locations 

Method 

Mults. 

Adda 


Trig. 

Tiig.*i 

Lat-Lon 


dn^-t- 127 m. 

8n,'t’81 m, 

S m, 

14 m, 

3 m. 

PROPOSED 

24 Op 106 nip 

13 Op + 53 nip 

Hp ♦ 6 nip 

2np+ 14mp 


Cartesian 


9n,-t- 138 m, 

12n, + 99m, 

4 m, 

14 m. 

m. 

PROPOSED 

19 Op + 101 nip 


3mp 


2nip 


Table S 

Comparison of Methods for an Actual Situation 


Coordinate System Desired 
for Locations 

Method 

Mults. 

Adds 

>r 

Trig. 

Trig.*‘ 

Lat-Lon 

SPLINE 

3946 

2788 

140 

392 

84 

PROPOSED 

2396 

1289 

103 

210 

190 

Cartesian 

SPUNE 

4449 

3552 

112 

392 

28 

PROPOSED 

1931 

1196 

6 

210 

4 
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